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Q ■ Abstract 

In the fragmentation of a heavy quark into a heavy meson whose hght degrees 
of freedom have angular momentum 3/2, all the helicity probabilities are completely 
fT^ ■ determined in the heavy quark limit up to a single probability w-^/2- We point out 

^^ . that this probability depends on the longitudinal momentum fraction z of the meson 

0> ' and on its transverse momentum pj_ relative to the jet axis. We calculate W3/2 as a 






X 



function of scaling variables corresponding to z and pj_ for the heavy quark limit of the 



Q^' perturbative QCD fragmentation functions for b quark to fragment into (be) mesons. 

rS^ • In this model, the light degrees of freedom prefer to have their angular momentum 



aligned transverse to, rather than along, the jet axis. Implications for the production 
of excited heavy mesons, like D** and B**, are discussed. 



The discovery of the powerful heavy quark spin-flavor symmetry [Q in Quantum 
Chromodynamics (QCD) and the development of the Heavy Quark Effective Theory (HQET) 
[^ have greatly improved our theoretical understanding of hadrons containing a single heavy 
quark Q. The crucial idea in HQET is that both the heavy quark mass and spin decouple 
from the strong interaction dynamics between the heavy quark and the light degrees of 
freedom in the limit of infinitely heavy quark mass. The decouplings of the mass and spin 
occur because in the limit rriQ -^ oo the only variable to describe a free propagating heavy 
quark is its 4-velocity v. In the presence of gauge fields, the propagation of the heavy quark 
is simply described by the Wilson line Pexp i J A^v^dt, which contains no information about 
the heavy quark flavor nor its spin. Furthermore, the leading operator - the chromo-magnetic 
dipole moment that couples heavy quark spin to the gluon field is inversely proportional to 
the heavy quark mass. To the first approximation in the heavy quark mass expansion, a 
heavy quark can be treated as a static color source for the remaining light degrees of freedom 
that make up the physical heavy-light hadrons observed in Nature. 

Falk and Peskin have recently pointed out that heavy quark spin symmetry pro- 
vides strong constraints on the helicity probabilities of the heavy mesons produced by the 
fragmentation/hadronization of a heavy quark. For heavy mesons whose light degrees of 
freedom have angular momentum 1/2, the helicity probabilities are completely determined 
in the heavy quark limit. If the light degrees of freedom have angular momentum 3/2, the 
probabilities are determined up to a single parameter W73/2. In this paper, we point out that 
the probability W3/2 depends on the longitudinal momentum fraction z of the heavy meson 
relative to the heavy quark jet and its transverse momentum p_|_ relative to the jet axis. 
We calculate W3/2 as a function of the scaling variables corresponding to z and p±_ for the 
nib —>■ 00 limit of the perturbative QCD (PQCD) fragmentation functions for a b quark to 
fragment into (be) mesons. These fragmentation functions can be used as a model for the 
fragmentation of a heavy quark into heavy-light mesons. The implications for the production 
of excited heavy-light mesons, like the D** and B**, are discussed. 

In the heavy quark limit it is convenient to label hadronic states by the eigenvalues 
j and ji of the total angular momentum J of the hadron and the angular momentum J/ of 
the light degrees of freedom respectively M. In general, the spectrum of hadrons containing 



a single heavy quark Q has, for each j;, a degenerate doublet (j_,j+) with total angular 
momentum j± = ji ± 1/2. (For the exceptional case of a baryon with ji = 0, there is only 
a singlet with total angular momentum 1/2.) For the heavy-light (Qq) mesons, one can 

write Ji = Sq + L where 5*^ is the spin of the valence light antiquark q and L is the orbital 
angular momentum with integer eigenvalue L. For S-wave heavy-light mesons, L = and 
ji^ = I , we have the familiar degenerate doublet (j_^,j+^) = (0~, 1~) consisting of a 
pseudoscalar and a vector meson. For P-wave heavy-light mesons, L = 1 and the ji^ of 
the light degrees of freedom can be either | or | . In this case, we have two distinct 



doublets {j_^,j+^) = (0+, 1+ ) and (1+, 2+) for ji^ = ^ and | respectively. While the 
0^ and 2~^ states are identified as the ^Pq and ^P2 states constructed in the LS coupling 
scheme respectively, the 1^ and 1^ states are linear combinations of the ^Pi and ^Pi states: 

|1+') = yiTsl^Pi) + v^pPi) and |1+) = -v/273|iPi) + ^fPr) . 

Falk and Peskin [Q showed that heavy quark spin symmetry combined with the parity 
invariance of QCD interactions can impose very useful relations among the probabilities for a 
heavy quark with a given helicity to fragment /hadronize into the various helicity states of the 
heavy-light hadrons within the same doublet. For definiteness, we assume the heavy quark 
Q is purely left-handed in what follows. For the case of ji = |, parity invariance implies 
that the two fragmentation probabilities for the heavy quark Q to hadronize by combining 
with light degrees of freedom with helicity rrii = —1/2 or +1/2 must be the same. Although 
heavy quark symmetry allows coherent superposition of the two zero helicities states in 
the corresponding (0, 1) doublet at the very early stage of the fragmentation process, the 
relative population of the four helicities states in the doublet are completely determined by 
symmetry. The resulting table of fragmentation probabilities is |^ 

i i \ 

2 4 ^ \ 

J' 

where the helicity h runs through the values —1, 0, +1 across the table. However, for ji = |, 
parity invariance implies the following table of fragmentation probabilities for a heavy quark 




to fragment into various helicity states of the light degrees of freedom [Q 

P3/2{mi) = I^-Wy2, -(1 - W3/2), -(1 - ^3/2), -W3/2J , (2) 

where the hehcity mi of the hght degrees of freedom runs through the values — |, — |, +|, 
+ 1 across the table. The Falk-Peskin parameter, W3/2, is the conditional probability for 
the heavy quark to fragment into a meson system whose light degrees of freedom are in 
the maximum helicity states mi = ±3/2. It can take values between and 1. At the 
very early stage in the fragmentation process, heavy quark symmetry allows coherent linear 
superposition of the various components in the (1, 2) doublet which have the same helicity 
h = m,i — ^. However, at a later time determined by the mass splittings within the doublet, 
the various helicity components of the (1,2) doublet will propagate incoherently, resulting 
in the following table of fragmentation probabilities 

]w3/2 §(1-^3/2) \{l-W3/2) 1^3/2 

(3) 




l{l-W3/2) \{l-W3/2) 1^3/2 

with the helicity h runs through the values —2, —1, 0, +1, +2 across the table. Therefore, 
all the relative fragmentation probabilities for the 8 helicity states of a (1, 2) doublet are 
determined, up to the Falk-Peskin parameter 103/2- This non-perturbative parameter can 
not be determined without a dynamical calculation. 

Consider the leading order Feynman diagram shown in Fig.l for a heavy b quark to 
fragment into a heavy-heavy (be) meson by creating a cc pair from the vacuum. In the limit 
of mfy/rric — > 00, this can be taken as a model for heavy quark fragmentation into heavy-light 
mesons. Following our previous works in g, |^, (see also 0), the fragmentation function 
D{z) can be expressed as D{z) = J^^^^-^ dsD{z, s), with so{z) = M'^/z + r^M^/(l — z) and 

D{z,s) = {16TT'^)~^\im.g^/mi-*oo{\-M.\'^/\-M.o\'^)- M. is the amplitude for a high energy source 
(symbolically denoted by F in Fig.l) to create b*{q) -^ {bc){p) + c{p') with total 4- momentum 
q = p + p' ., and M.q denotes the same source creating an on-shell heavy b quark with the 
same 3-momentum q. In the nonrelativistic approximation, M = m,f, + mic', and we define 
r = rric/M. s = g^ is the virtuality of the fragmenting b quark. In a frame where the virtual 



b* quark has 4-monientum q = (go, 0,0,^3), the longitudinal momentum fraction of (be) is 
^ = (Po + 'P^)l{% + 53) aiid its transverse momentum is pj_ = (^1,^2)- For the b quark 
propagator and bbg vertex that entered in Fig.l, we will apply the techniques of HQET |^. 
It is straightforward to write down the following matrix elements M.{}Pi^ h) and Ai(^Pi, h) 
for b -^ (be), with the (be) bound state in definite helicity h of the ^Pi and ^Pi states 
respectively: 



MCP.,k) = u,,^c,RXo)^^^l^^^,Mm-,ry' (ii^j r . (4) 

with the vertex V^ given by 

V- = ^rMq^ fi _ !i:A ^^ 2rM^^n" ^ + ^il + ^ i]^'^ ■ (5) 

y n ■ k J [n ■ ky \ n ■ k J 

and 



with the vertex tensor V^^ = V^j^, where V^ can be obtained from V"' in (|^) with fi ^ — fi 
and 7" -^ —7". In (^ and (|^), i and j are the fundamental color indices, Nc is the number 
of color, Cf = {N^ — l)/{2Nc), and -R'(O) and e{h) are the derivative of the radial wave 
function and the helicity wave function for the bound state respectively. We have picked the 
axial gauge associated with the vector n^ = (1,— p7|p|). In this gauge, the short- distance 
process that creates the energetic heavy quark b* (symbolically denoted by F in Fig.l) and 
the subsequent fragmentation of the heavy quark become manifestly factorized. In the 
above equations, we have used p = Mv and q = mf,v + k, where v and k are the 4-velocity 
and residual momentum of the b quark respectively. Some useful kinematical relations are 
fc2 = 2rv ■ k = r{s — (1 — rYM'^)/M. The tree-level matrix element squared |A1oP is given 
by A^c(^/-2)Tr[rr(l + /;)]. In writing down the above matrix elements, we have assumed 
the P-wave bound state is a color-singlet. It has been pointed out recently by Bodwin- 
Braaten-Lepage that, beyond leading orders in the calculations of the production and 
decay rates of P-wave quarkonium [^, |T0[ and (6c) mesons []^, there is also a color-octet 



S-wave mechanism needed to be taken into account in order to avoid infrared divergencies 
that spoil factorization. We will also include the color-octet S-wave contributions in what 
follows. 

Heavy quark fragmentation functions have nontrivial limits as rriQ —>■ oo (or r ^ 0) if 
they are expressed in terms of the scaling variable y = iXj z — l + r)/r, first used by Jaffe and 
Randall |jll|, and the rescaled transverse momentum t = \p±\/{rM). We therefore define 



the fragmentation functions D{y,t), D{t), and D{y) according to the following changes of 
variables: 

1 rl POD 

dzD{z) = dz dsD{z,s) , (7) 

Jo Jso{z) 

("OO 

dy dtD{y,t), (8) 

1 JO 

oo 





dtD{t), (9) 

1 

dyDiy). (10) 

The relation among s, t, and z is given by s = M^[(l + r'^t^)/ z + r^(l + t^)/(l — z)]. Fol- 
lowing the same procedures as in Refs.[^ ^ 0], we square the matrix elements M.{}Pi,h) 
and A^('^Pi, h), calculate the interference between these two amplitudes, and project out the 
transversely and longitudinally polarized l"*" states, one can deduce the generalized polar- 
ized fragmentation functions Dxiy^t) = D^\y,t) + 2D^^\y,t) and DL{y,t) = D]^\y,t) + 
D^^\y,t), depending on both variables y and t. In the limit r — > 0, the color-singlet pieces 
DT^iVit) are given by 

D^T\y.t) = ^ yi^^^^y!^, {{y - 4) V + 1/^(40 - 32y + 26y' - V + v'^ 

+ 3y\l6 -8y + 2y^ + 3y^)t^ + 3y\8 - Ay + 3y^)t'' + (1 + y)¥} , (11) 

D'^{y,t) = ^4£f^((i/-4)V + y^(4 + 40y-10y^-V + yV 



3 y\t^ + y^f 
+ 3y\A + I3y - Ay^)t^ + 3y\h + 2y)t^ + (1 + yft^] , (12) 



where Ci = 2a'^\R'(0)\'^/(3TcNc'ml). The color-octet piece D^^\y,t) can be extracted from 
Ref.|T2| by taking the hmit r ^ 0, 



D^'\y,t) = 12Cs 



{y - m 



y^{t^ + y 



2M 



%2 + y{y + 4)t2 + t' 



(13) 



with Cg = 3a;^ifg/(256A^c"^c) where H'^ is a nonperturbative parameter associated with the 
color-octet mechanism for P-wave production [0, ^ |TU]. We will treat the overall constants 
Ci and Cg as free parameters in our approach. One notices that in terms of the variables y 
and t, the leading order results of the fragmentation functions DT,L{y,t) scale, i.e. they do 
not depend on the heavy quark mass. 

This motivates the introduction of the generalized Falk-Peskin parameter w^i2{y,t) 
that is a function of the two scaling variables y and t. The generalized polarized fragmenta- 
tion functions D{y, t) for the 8 helicity components of a (1, 2) doublet satisfy a similar table 
like that given in (^) with W3/2 replaced by W3/2{y,t). From the probability table for Pi{h) 
in (^, we see that Dt is proportional to (1 — W3/2)/8 + 81^3/2/8, while Dl is proportional 
to (1 — ^3/2)74. Thus W3/2{y,t) can be defined in terms of the fragmentation functions for 
the spin-1 state by 

DT{y,t)-lDLiy,t) 



W3/2{y,t) - „ . +\ , n f +\ ■ 
DT{y,t)+DL{y,t) 

Simple analytic expression for W3/2{y,t) can be obtained from ([TT|)-(pr3D as 



(14) 



W3/2{y,t) 



n{y,t) + €h{y,t) 
d{y,t) + 2eh{y,t) 



(15) 



withn{y,t) = 6y\y-l)H^[Ay^ + y{y + A)t^ + t''],d{y,t) = 2{t'' + y^)[y\y-4f + y^{24 + y 
Ay^ + y^)e + y\n -2y + 2y^)t^ + (1 + y^t"^] , h{y, t) = 9y\t^ + 1/2)2 [4^2 ^ y^y ^ 4)^2 ^ ^4 
and e = Cg/Ci. 

We can define probabilities u'3/2(t) and 1^3/2 (y) that depend only on a single scaling 
variable by an expression analogous to (|T^), except with Dx^Liy^t) replaced by Dx^iii) and 
DT,L{y) respectively. Integrating the fragmentation functions in (pT])-([T3|) over y and forming 



a similar ratio like ([I^), we obtain wz/2{t): 



9 ( rii + n2 arctan(t) + ^3 log(l + t^) 

^3/2(t) = - 



80 \ (ii + (^2 arctan(t) + d-^ log(l + t^) 



(16) 



with rii = t 630 - 5 (521 - 240 e) t^ + (231 - 2440 e) t^ 



n2 



126- 15(49- 16e)t2 + 



8(20-99e)tH(5 + 72e)t6 , rig = -320t 4-2 (2-3e) t2_3et4 ^ ^^^ ^^ ^^ 105-2(406- 



135 e) t^ + (79 _ 549 e) t^ 



d. 



35-5 (71 - 18 e) t^ + 3 (31 - 99 e) t^ + 3 (1 + 9 e) t^ 



d. 



-72 1 



4-6(l-e)t2 -3et^ 



The curve of 1^3/2 (t) versus t is plotted in Fig. 2 for the two 



cases of e = (color-singlet dominance) and e = 00 (color-octet dominance). 



Similarly, by integrating the fragmentation functions in (|Tl])-(|T3|) over t and forming 
the ratio like (0), we obtain w^i2{y)'- 



wz/2{y) 



1 {y - 1)2(12 + 81/ + 5i/2) + 15 e 1/2(8 + Ay + 3y^) 
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(8 + 4|/2 + y4) + sey^{8 + 4:y + 3y^ 



(17) 



The curve of W3/2{y) versus y is plotted in Fig.3 for the two cases of e = and 00. 

The original Falk-Peskin parameter ^3/2 is given by (0), with the numerator and de- 
nominator integrated over both y and t. The total fragmentation probabilities are J^ dy /q°° dtDxiy, t) 
(86/315)Ci + (16/5)^8 and J^ dy J^ dtDL{y,t) = (17/63)Ci + (8/5)^8. These imply 



^^3/2 



1 /29 + 504e 
6 U9 + 168e 



We note that setting e = in ([TsD gives U73/2 — ^ 29/114, a result first derived by Chen 



and Wise [13]. While the Falk-Peskin parameter 1^3/2 is independent of QCD evolution, its 
generalizations w^/2{y,t), Ws/2{t), and ^3/2(1/) given above are evaluated at a subtraction 
point around the heavy quark mass scale. Their values at a higher scale can be determined 
by the usual Altarelli-Parisi evolution of the fragmentation functions. 

The Falk-Peskin parameter 1^3/2 determines all the helicity probabilities for the frag- 
mentation of a heavy quark into a heavy meson whose light degrees of freedom have angular 



momentum 3/2. In this Letter, we have pointed out that W3/2 can have nontrivial depen- 
dence on the scahng variables y and t corresponding to the longitudinal momentum fraction 
z and the transverse momentum p±_ of the meson relative to the heavy quark jet. We have 
calculated W3/2 as a function of y and t for the m,6 — > 00 limit of the perturbative QCD 
fragmentation functions for the production of P-wave (6c) mesons. We find that the proba- 
bility U73/2 has significant dependence on y and t, varying over the range from to 1/2. The 
PQCD fragmentation functions can be used as a model for the fragmentation of a heavy 
quark into heavy-light mesons, and applied to the fragmentation processes c — > D** and 
h -^ B**. The prediction ws/2 < 1/2 of this model implies that light degrees of freedom 
with helicity mi = ±1/2 always have a larger population than the maximum helicity states 
of mi = ±3/2. This prediction supports the speculation of Falk and Peskin that the 
angular momentum of the light degrees of freedom with ji = 3/2 prefers to align transverse 
to, rather than along, the fragmentation axis. This spin alignment of the heavy quark can 
be detected by measuring the anisotropics of the decay products from the hadronic transi- 
tions between the two doublets (l+,2+) and (0^, 1^) Q. The PQCD fragmentation model 
predicts that these anisotropics vary significantly with y and t. In the case of charm quark 
fragmentation into D**, an upper bound of 1^3/2 < 0.25 at the 90% confidence level has 
been deduced from the existing experimental data B. After integrating over y and t, the 



prediction (|T8|) of the PQCD fragmentation model with e = is W3/2 ~ 0.25, which suggests 
that present experiments may be close to observing a nonzero value for ^3/2- String models 
of fragmentation tend to give significantly larger values of W3/2 and may already have been 
ruled out [Q. We conclude that future measurements of the Falk- Peskin probability W3/2 for 
the charm and bottom systems, including the dependence of W3/2 on the scaling variables y 
and t, can provide valuable insights into the dynamics of heavy quark fragmentation. 

I am grateful to Eric Braaten for valuable conversations and Mark Wise for useful 
communication. This work was supported in part by the U.S. Department of Energy, Division 
of High Energy Physics under Grants DE-FG03-91ER40674. 
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Figure Captions 

1. Leading order diagram for the heavy b quark to fragment into a P-wave (be) bound 
state as a model for heavy quark fragmentation into excited heavy-light meson. The 

9 



outgoing momenta are (1 — r)p + 6, rp — 6, and p' for the b, c, and c, respectively. 6 is 
the relative momentum of the b and c. 

2- 'W3/2{t) versus the rescaled transverse momentum t at the heavy quark mass scale. 

3- u]3/2{y) versus the scahng variable y at the heavy quark mass scale. 
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